Dynamics of nematic elastomers 
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We study the low-frequency, long-wavelength dynamics of soft and semi-soft nematic elastomers 
using two different but related dynamic theories. Our first formulation describes the pure hydrody- 
namic behavior of nematic elastomers in which the nematic director has relaxed to its equilibrium 
value in the presence of strain. We find that the sound-modes structure for soft elastomers is 
identical to that of columnar liquid crystals. Our second formulation generalizes the derivation of 
the equations of nematohydrodynamics by the Harvard group to nematic elastomers. It treats the 
director explicitly and describes slow modes beyond the hydrodynamic limit. 

PACS numbers: 83.80.Va, 61.41. -|-e, 62.30.-|-d 



I. INTRODUCTION 

Nematic elastomers (NEs) Q, |^ |M IS are rubbery mate- 
rials with the macroscopic symmetry properties of ne- 
matic liquid crystals i3j. In addition to the elas- 
tic degrees of freedom of ordinary rubber, nematic elas- 
tomers possess the internal, orientational degree of free- 
dom of liquid crystals. Since NEs are amorphous solids 
rather than fluids, their mechanical properties differ sig- 
nificantly from those of standard nematic liquid crystals. 
The interplay between elastic and orientational degrees 
of freedom is responsible for several fascinating proper- 
ties of NEs. For example, temperature change or illu- 
mination can change the orientational order and cause 
the elastomer to extend or contract by several hundred 
percent l a la B . Nematic elastomers display a soft elas- 
ticity ITllll^ characterized by vanishing shear 
stresses for a range of longitudinal strains applied per- 
pendicular to the nematic direction. They exhibit an 
anomalous elasticity in which certain bending and shear 
moduli are length-scale dependent and vanish or diverge 
at long length-scales 

In addition to unusual static properties, NEs have an 
intriguing dynamic mechanical behavior. Earl y rh eology 
experiments on liquid crystalline elastomers [l7j found 
no nematic effects in NEs. Later, however, several ex- 
periments [ISj, (ij, I23, 01 observed a genuinely uncon- 
ventional response to oscillatory shear. It turns out that 
an internal relaxation of the nematic director leads to a 
dynamic mechanical softening of NEs. This behavior has 
been christened dynamic soft elasticity. It makes NEs in- 
teresting for device applications in areas such as mechan- 
ical vibration damping (exploiting the fact that the 
mechanical loss is record-high over a wide range of tem- 
peratures and frequencies) or acoustics where NEs open 
the possibility of acoustic polarization (using the fact 
that only particular soft shears are strongly attenuated) 
analogous to the optical polarization in birefringent me- 
dia. 

The theoretical investigation of the dynamic- 
mechanical properties of NEs was pioneere d by Terentjev, 
Warner and coworkers (TW) 0, 111 IH IM Hi. Ref- 
erences Q and ll^l present a detailed derivation of 



equations governing the long-wavelength low-frequency 
dynamics of NEs and derive their associated mode 
structure, which includes non-hydrodynamic, rapidly 
decaying modes. This derivation, which assumes a 
single relaxation time for the director, is based on the 
Lagrangian approach 2^ to the dynamics of continuum 
systems in which non-dissipative forces are calculated as 
derivatives of an energy functional and dissipative forces 
as derivatives of a a Rayleigh dissipation function. The 
TW work focuses primarily on rheological response in 
both soft and semi-soft NEs at zero wavenumber, and 
it neglects contributions to dynamical equations arising 
from the Frank free energy for director distortions, which 
are higher order in wavenumber than those arising from 
network elasticity. As a result, as we shall see, the mode 
structure derived in Refs. 0] and misses diffusive 
modes along symmetry directions in soft NEs. When 
contributions from the Frank free energy are included 
and those from the elastic network are excluded, the 
TW approach reproduces the original Leslie-Ericksen 
equations .26 , 27] of nematodynamics, excluding the 
rotational inertial term that is usually discarded [23 or 
is missing entirely in alternative derivations Q of these 
equations. 

In this paper, we will present alternative approaches 
to deriving the equations governing the dynamics of NEs. 
First we will derive the exclusively hydrodynamical equa- 
tions for those variables whose characteristic frequencies 
vanish with wavenumber. These equations, like those for 
the hydrodynamics of smectic 0,129] and columnar liq- 
uid crystals, exclude non-hydrodynamic director modes. 
Interestingly, the modes for soft NEs predicted by these 
equations are identical in form to those of columnar liquid 
crystals W\ with three pairs of sound modes and diffusive 
modes along symmetry directions where sound velocities 
vanish. One pair of sound modes is predominantly lon- 
gitudinal with non-vanishing velocity at all angles. A 
second pair is predominantly transverse with a velocity 
that vanishes for wavevector either parallel or perpendic- 
ular to the uniaxial symmetry axis, and a third pair is 
completely transverse with a velocity that vanishes for 
wavevector parallel to the symmetry axis but is nonzero 
otherwise. We then derive the phenomenological equa- 
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tions for the slow dynamics of all displacement and direc- 
tor variables in NEs using the Poisson-Bracket formalism 
[iSI^ for obtaining the dynamics of coarse-grained vari- 
ables, applied so successfully to the study of dynamical 
critical phenomena 32J, for deriving phenomcnological 
equations for any set of coarse-grained variables whose 
dynamics is slow on a time scale set by microscopic colli- 
sion times. Like TW, we assume a single relaxation time 
for the director. We do, however, discuss how this con- 
straint can be relaxed. When applied to fluid nematic 
liquid crystals, the Poisson-bracket formalism is equiva- 
lent to that used by the "Harvard group" [s^, in its 
derivation of hydrodynamics of nematics. Our dynamical 
equations for NEs reduce to the equations of nematohy- 
drodynamics derived by the Harvard group when elastic 
rigidities vanish and to the purely hydrodynamics equa- 
tions for NEs when fast modes are removed. When con- 
tributions from the Frank free energy are ignored, our 
equations are identical to those derived by TW. We will 
not discuss random stresses or inhomogeneities in this 
paper, though they contribute static components in light 
scattering experiments 35| that may obscure the obser- 
vation of the NE modes we discuss here. 

The outline of our paper is as follows. In Sec. ITU we 
briefly review the general Poisson-bracket formalism for 
obtaining coarse-grained dynamics. Section IHII focuses 
on the pure hydrodynamics of NEs. We set up equa- 
tions of motion for the momentum density and the elastic 
displacement. We derive the appropriate elastic energy 
entering these equations by integrating out the director 
degrees of freedom. Then we compare the hydrodynamic 
equations to those of conventional uniaxial solids and 
columnar liquid crystals. After extracting the sound ve- 
locities of the modes we finally determine the full mode 
structure in the incompressible limit. Section lTvl features 
our formulation that explicitly accounts for the dynam- 
ics of the director. We set up equations of motion for 
the momentum density, the elastic displacement and the 
director. These are then compared to the equations of 
motion for uniaxial solids and to the equations of motion 
derived by TW. We determine the mode structure. Fi- 
nally, we compare the so obtained modes to the modes 
obtained by TW. A brief summary is given in Sec. 



II. COARSE-GRAINED DYNAMICS 

Stochastic dynamical equations for coarse-grained 
fields [31I Is^ l can be obtained by combining the Poisson- 
Bracket formalisms of Classical mechanics |23, which 
guarantees the correct reactive couplings between fields 
with opposite signs under time reversal, and the Langevin 
[s^ approach to stochastic dynamics, which provides a 
description of dissipative processes and noise forces. Let 
<l)^;(x,t), ^ = 1,2,... be a set of coarse-grained fields 
whose statistical mechanics is described by a coarse- 
grained Hamiltonian TL. The dynamical equations for 



<I>^(x,t) are first-order differential equations in time: 



5n 



+ Cm(x, t). 



(2.1) 



Here and in the following the Einstein summation con- 
vention is understood. The first term on the right 
hand side is a non-dissipative velocity, also known as 
the reactive term, that contains the Poisson bracket 
{<i>^(x, t), $,y(x', t')} of the coarse grained fields W^. The 
reactive term couples to ST-L/S^y only if <i>^ and 
have opposite signs under time reversal (when external 
magnetic fields are zero). The second term on the right 
hand side is a dissipative term. F^^j/ are the components 
of the so-called dissipative tensor. This tensor couples 
<i>^ to SH/S^i, only if <I>^ and have the same sign 
under time reversal. If the noise term is present, 
Eq. H2.1|l represents a stochastic or Langevin equation. 
As such E g. (12.111 may be used to set up a dynamic func- 
tional [sM I39I l40l| to study the effects of nonhnearities 
and fluctuations via dynamical field theory. In this pa- 
per we are not interested in these effects. Hence, we focus 
on linearized hydrodynamic equations and pay no further 
attention to noise. 



III. PURE HYDRODYNAMICS 

Hydrodynamics describes the dynamics of those de- 
grees of freedom whose characteristic frequencies uj van- 
ish as wavenumber tends to zero. In other words 
hydrodynamics focuses exclusively on the leading low- 
frequency, long-wavelength behavior. All other degrees 
of freedom, even though they might be slow, are strictly 
speaking not hydrodynamic ones. In this section we will 
derive the hydrodynamic equations of NEs. These equa- 
tions apply for frequencies lu such that ujt <ti 1 , where r is 
the longest non-hydrodynamic decay time in the system. 
As we will show in the next section, the characteristic 
time for director decay is in fact very slow with r ^ 10^^ 
sec. 0, so that the regime of applicability of hydro- 
dynamics is quite small for current NEs. It is imaginable, 
however, that other systems will be found with shorter 
decay times. 

There are two general classes of hydrodynamic vari- 
ables: conserved variables and broken symmetry vari- 
ables. A single-component NE has the same set of con- 
served variables as a fluid: energy density e , mass den- 
sity p, and momentum density g. It also has the same 
broken symmetry variables as a crystalline solid, namely 
three displacement variables, though strictly speaking, 
these variables in an elastomer are not associated with 
a macroscopic broken symmetry because in the absence 
of orientational order elastomers have the same macro- 
scopic rotational and translational symmetry as a fluid. 
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Since it is rotational symmetry that distinguishes a ne- 
matic elastomer from an isotropic one, it could be argued 
that the nematic director should be a hydrodynamic vari- 
able, but as in smectic and columnar liquid crystals, the 
director degrees of freedom decay in microscopic times 
to their preferred configuration in the presence of strain 
and are thus not hydrodynamic variables. Elastomers 
differ from equilibrium crystals in at least two impor- 
tant ways: The first, alluded to above, is that they are 
not periodic and thus do not have mass-density wave or- 
der parameters whose phases act as broken symmetry 
hydrodynamic variables; rather Lagrangian displacement 
variables u take their place. The second is that an elas- 
tomer is permanently crosslinked: it is a classical rubber 
in which changes 5p in mass density are locked to changes 
in volume such that 5p/ p = — V-u and in which perme- 
ation in which there is translation of a mass-density- wave 
without mass motion is prohibited. Thus, mass density 
is not an independent hydrodynamic variable, and we are 
left with a total of 5-1-3— 1 = 7 independent variables and 
7 associated hydrodynamic modes. These modes are heat 
diffusion and, depending on direction, either six propa- 
gating sound modes, four propagating sound modes and 
two diffusive displacement/ velocity modes, or two longi- 
tudinal sound modes and diffusive displacement/ velocity 
modes. In what follows, we will consider only isothermal 
processes so that heat heat diffusion can be ignored. 



A. Elastic energy 

To derive the hydrodynamical equations for NEs, we 
first need the appropriate elastic free energy. The elas- 
tic constant measuring the energy of strains in planes 
containing the anisotropy axis [the shear modulus C5 to 
be defined in Eq. (|3.2I) ] vanishes as a result of the broken 
symmetry brought about by the establishment of nematic 
order. Thus a good starting point for this free energy is 
that of a uniaxial solid with this elastic constant sim- 
ply set to zero. This leaves certain soft directions in 
which distortions cost zero energy, and as in smectic and 
columnar liquid crystals 5] , curvature- like terms that are 
quadratic in second-order spatial derivatives have to be 
added to ensure stability. To make contact with dynam- 
ical equations involving the director to be presented in 
Sec. IIVI rather than simply adding these terms, we find 
it useful to derive them from the free energy of a nematic 
elastomer expressed in terms of both the strain and the 
directors. We will restrict our attention to harmonic dis- 
tortions. 

Elastomers are permanently crosslinked pieces of rub- 
ber whose static elasticity is most easily described in La- 
grangian coordinates in which x labels a mass point in the 
unstretched (reference) material and R(x) — x-l-u(x), 
where u(x) is the displacement variable, labels the posi- 
tion of the mass point x in the stretched (target) mate- 
rial. We will use Lagrangian coordinates throughout this 
paper. We will, however on occasion make reference to 



Eulerian coordinates in which r = R(x) specifies a posi- 
tion in space and u(r) the displacement variable at that 
position. 

The elastic energy of a nematic elastomer can be di- 
vided into three parts: 



(3.1) 



where Hu is the usual elastic energy of a uniaxial solid. 
Tin is the Frank free energy of a nematic, and 7iu,n is 
the energy of coupling between strain and director dis- 
tortions. 

Choosing the coordinate system so that the z-direction 
coincides with the uniaxial direction, we have to har- 
monic order lial 



C C 



(3.2) 



Here, linearized components of the Lagrange 

strain tensor u: 



1 



- {d,uj + djUi) . 



(3.3) 



We will use the convention that indices from the be- 
ginning of the alphabet, {a, 6}, assume the values 1 
and 2 whereas indices from the middle of the alphabet, 
fc, Z}, run from 1 to 3. Note that, com par ed to our 
work on the anomalous elasticity of NEs [l^ll5l| the elas- 
tic constants in Eq. (|3.2(l have a somewhat different def- 
inition that is geared towards taking the incompressible 
limit. Here we have arranged things so that the terms 
featuring C2 and C3 involve the trace of the strain ten- 
sor. In the incompressible limit that we will eventually 
take one has uu = so that C2 and C3 drop out. 

Expanded to harmonic order in the deviation 6n = 
n — no from the uniform equilibrium state Uq — e^, the 
Frank energy reads 



'Hsn — 



^ {dMa) 



(3.4) 



where Sab = —£ba is the two-dimensional Levi-Civita 
symbol. For notational simplicity in the remainder of 
this paper we will replace 6n by n with the understanding 
that it has only two components Ua- With this notation, 
Hsn can be expressed as 

nsn = -^ J (fxnaMabiV)nb, (3.5) 

where 

Mab = {Ki - K2)dadb + {K2dl + K:idl)6ab (3.6) 
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and where 9^ = dada with the Einstein convention un- 
derstood. 

The couphng energy, finally, can be written in the form 



Wu,„ = J d^xl^^Ql+D2U,aQa 



(3.7) 



where once more terms beyond harmonic order have been 
neglected and where 



1 



[dzUa - daUz) = 5na - Via 



(3.8) 



As explained above the director is not a genuine hy- 
drodynamic variable and hence it should be integrated 
out of the elastic energy as long as we focus on the hy- 
drodynamic limit. We do so by minimizing 7i over Ua to 
find 



- 6 ^2 

^1 



6 ^2 
^1 



(3.9) 



Inserting this equation into Ti. and retaining only the 
dominant terms in gradients, we obtain 

/r C C 



where 



ci 



^3 



C5 

i 
i 



1 - 



Dl 
2Di ' 

Dl 
D2 



Ki 



(3.10) 

(3.11a) 
(3.11b) 
(3.11c) 



The superscript R indicates that C5, Ki and have 
been renormalized by director fluctuations. Note that 
the elastic constants Kf- and are not the same as 
the Frank splay and bend constants Ki and K^. This 
is in contrast to smectic-A and columnar liquid crystals 
in which the coefficients of quartic gradient terms in the 
effective elastic energy arising from the relaxation of di- 
rector modes are identical to the Frank elastic constants. 
The modulus D2 can have either sign, and there are no 
thermodynamic constraints preventing cither 1 -|- D2 / Di 
or 1 — D2/D1 from being zero. Thus, either or 
in Eq. H3.11|l could be zero. In this case, contributions to 
(9^Mz)2 or {dIuaY in '^u-ci arising from network elas- 
ticity would have to be added to ensure stability. Except 
for the exceptional cases when D2 — ±Z?i, those contri- 
butions are generally smaller than the ones represented 
in Eq. H3.11|l . and we will ignore them. If C§ vanishes, 
then the dependence of the elastic energy on u^z drops 
out. This is the famous soft elasticity. If the condition 



for softness (|3.11a|l is not strictly fulfilled and there is 
a small but finite remnant then small shears in the 
planes containing the anisotropy axis do cost a small en- 
ergy. This non-ideal behavior is known as semi-softness. 

The energy (|3.10l) combines attributes of both smec- 
tic and columnar liquid crystals. The vertical displace- 
ment Uz is analogous to the displacement variable u of 
a smectic liquid crystal. It is soft for distortions in the 
_L-direction, and a bending term iff- (9^ Uz) 2/2 is needed 
to stabilize it. The in-plane displacements Ua are analo- 
gous to those of columnar liquid crystal, which are soft 
for distortions in the z-direction, and a bending term 
K^\d1uaf' 12 is need to stabilize them. 



B. Hydrodynamic equations 

We are now in a position to write down the full hy- 
drodynamic equations for nematic elastomers. For sim- 
plicity, we will restrict our attention to isothermal pro- 
cesses so that we can ignore temperature diffusion. This 
leaves us with six hydrodynamical variables, the momen- 
tum density gi(x) and the displacement Mi(x). These 
are independent variables at the reference point x that 
satisfy the continuum generalizations of the usual rela- 
tions for the momentum and displacement of a parti- 
cle: 9(7i(x)/35j(x') = (5y5(x — x'), dui{'x)/duj{'x') = 
6ijS{x — x'), and dgi{'x.) / duji'x') — 0. These relations 
yield a non- vanishing Poisson bracket between Ui (x) and 
5,(x') 



{wj(x),gj(x')} = %(5(x-x'). 

The g — g and u — u Poisson brackets are zero, 
these results, we obtain the equations of motion 



kin 



1 

-9i 
P 



5gi 

SHu-cl , ^ 

gi = h Vijkl OjdlVk ; 



(3.12) 
Using 

(3.13a) 
(3.13b) 



where v is the velocity field, and Tikin = / d^xgigi/ {2p) 
is the coarse-grained kinetic energy, rjiju is the viscosity 
tensor, which has five independent components. We can 
parametrize the stress tensor so that the entropy produc- 
tion from viscous stresses takes on the same form as the 
elastic energy Tiu, 



TS 



m .2 

2 



V2 UzzUu 



m .2 



(3.14) 



where the 77's are low frequency viscosities, rj^ is an ef- 
fective viscosity that is, as its counterpart C^, renormal- 
ized by director fluctuations. Equation H3.14|l contains all 
contributions to the entropy production equation in the 
truly hydrodynamic limit we are considering. Nonhydro- 
dynamic variables like the director have already relaxed 
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to their local equilibrium values. Thermodynamic stabil- 
ity requires rji > for i = 1, 3, 4, 77^^ > and 771773 > 77^. 

A few observations about Eqs. H3.13|l are in order. 
First they are identical in form to the equations for a 
conventional elastic material (without vacancy diffusion) . 
The distinction between such a material and a ncmatic 
elastomer appears only in the form of 7iu-ei- The ab- 
sence of any dissipative term proportional to —STi/Sui in 
Eq. (3.13a) reflects the tethered or crosslinked character 
of the elastomer. In non-crosslinked systems, this equa- 
tion would contain an additional dissipative term propor- 
tional to —6Tlu~ci/Sui describing permeation. Second, 
these equations are expressed in Lagrangian coordinates, 
and the derivatives u and g are time derivatives at con- 
stant value of the reference position x. In Eulerian co- 
ordinates, u and g become functions of points r — R(x) 
in space: U£;(r,t) — u(x(r),i) and similarly for g^;. The 
time derivative u — du/dt — dus/dt + v • Vu^ is the 
time derivative in the reference frame moving with the 
local fluid velocity. 

To obtain a more explicit form for the equations of 
motion, we note that Eq. (|3.13a|l allows us to replace gi in 
Eq. Ij3.13b|l with piii to produce the standard mechanical 
equation for a solid with dissipation, 

pila = (C2 + ri2dt)daUzz + (C3 + rizdt)daUii 

-K^dtua, (3.15a) 
pilz = (Ci + Tjidt + C2+ ri2dt)dzUzz 
+ (C2 + mdt + C3 + 'n3dt)dzUu 
+ {Cf + iiidt)dbUbz - K^dl^uz . (3.15b) 

Apart from the appearance of C|^, which is zero in a 
soft NE, rather than C5 and bending terms with quartic 
derivatives, these equations are identical to those of a 
uniaxial solid. 

Upon switching to frequency space one can see that 
the static renormalization of the shear modulus C5, 
Eq. (|3.11a(l . has a dynamical generalization with a fre- 
quency dependent C^{uj) = — ioJT]^, whose real part 
vanishes when = 0. This form for C^{uj) is only valid 
in the hydrodynamic limit with ujt ^ 1. As we will dis- 
cuss more fully in the next section, there are other terms 
in C^{uj) when this limit is not obeyed. 
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FIG. 1: Symmetry directions. The t-direction is perpendic- 
ular to the plane containing the equilibrium director no and 
the wavevector q of a generic excitation. The T-direction is 
perpendicular to q and the t-direction. 



The sound modes have frequencies a-'(q) = c{9) q where 
6 is the angle that q makes with the z-axis. The 
non-dissipative bending terms give rise to modes with 
Lo ^ along the symmetry direction. These modes, 
however, mix with dissipative ones and become over- 
damped diffusive modes with lu ^ —iq^. Thus, to ob- 
tain the true non-dissipative sound-mode structure, we 
can ignore the bending terms. Consequentially, the non- 
dissipative sound-mode structure is that of a uniaxial 
solid with C5 = 0, which it turns out is identical to that 
of columnar liquid crystals, though the input variables 
are slightly different (there is no in a columnar liquid 
crystal, but there is an independent density). The modes 
break up into displacements ut = Sab<laUb/q± perpendic- 
ular to both q and the z-axis and coupled displacements 
Ua/q± and in the plane that contains q (see 
Fig. nj. In Fourier space, the hydrodynamic equations 
are 

puP'Ut — C4 q\ut , (3.16a) 
puj'^u±_ = (C2 + C3) q^q^Uz + (C3 + 2C4) q'iu± , 

(3.16b) 

puj^Uz = (Ci + 2C2 + C3) qluz + (C2 + C3) q^qztli_ . 

(3.16c) 



C. Sound velocities 

To assess the mode structure of the equations of mo- 
tion H3.15|l . we begin with an analysis of propagating 
sound modes in the absence of dissipation. To keep ar- 
guments simple, we will restrict ourselves here to ideally 
soft NEs with = 0. When ^ 0, the bending terms 
can be neglected, and the dynamical equations and asso- 
ciated equations are identical to those of uniaxial solids. 
We will return to the case when is nonzero in Sec. lIVI 



Thus, there is a transverse sound mode with velocity 




(3.17a) 



that vanishes when 9 = and reaches a maximum at 
6 = tt/2. The sound velocities ci{d) and 02(0) of the 
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FIG. 2: Schematic polar plot (arbitrary units) of the sound 
velocities Ct (solid line), ci (short dashes), and C2 (long 
dashes) . 



other modes are coupled and satisfy 



cl cl = ^[(2^4 + C3)(Ci + 2C2 + Cs) 
P 



(3.17b) 



2 I 2 
"1 +C2 



= - [(2C4 + C3) sin^ e+{Ci+ 2C2 + C3) cos' d] 
P 

(3.17c) 



These equations are identical to the equations satisfied by 
the sound velocities in a columnar system ^ . One of the 
sound modes is purely longitudinal in the limit C3 — > 00; 
its velocity is non-vanishingfor all 9. The second mode 
is like that of a smectic-A with a sound velocity 

that vanishes at 6* = and 9 = tt/2. The third mode is 
a purely transverse sound mode whose velocity vanishes 
only at 6* = 0. Figure [3 plots the three sound velocities. 



D. Full incompressible mode structure 



Having found the general sound mode structure in the 
non-dissipative limit, we turn to a full analysis of modes 
in the incompressible limit. To discuss this limit, it is use- 
ful to decompose u(q) into a longitudinal part it/ along q 
and components Ut, perpendicular to both q and z, and 
ut perpendicular to q in the plane containing q and z as 
shown in Fig. ^ 

As was the case in the dissipationless limit, ut decou- 
ples from Ut and u/. In the incompressible limit ui van- 



ishes and we are left with 

2 , T^R„4: 



(Cl + 2C4) ^ + 4 + 4 



(3.18a) 
Ut 

Ut ■ 
(3.18b) 



These equations produce propagating modes with respec- 
tive frequencies 



ujt ± = ± 




(3.19a) 



P q 

2(771 + 2774) qlql + Tli [ql 
Apq^ 



(3.19b) 



as long as q is not along a symmetry direction in which 
the sound velocity is zero. When q± — 0, the t- and 
T-modes become diffusive with identical frequencies 



1 

Ap 



(3.20) 

As in conventional nematics Ki ^ ^ 10~^ dynes and 
p ^ 1 gram/cm'^. The viscosity should be larger than the 
0.01 poise characteristic of fluids. Thus, we can expect 
that 16pK^ ^ ivi')'^- In this limit the above modes 
become slow and a fast diffusive mode with frequencies 



Wt,s = LUt,s = -« 



— f 

V5 



■ V5 2 



(3.21a) 
(3.21b) 



When Qz = 0, the t-modes remain propagating sound 
modes, but the T-modes become diffusive with frequen- 
cies 



ujt,± 



1 

Tp 



iili±J~{rj^y + 16pK^ 



which in the limit 16K^p <C ivs')^ reduce to 



.2K^ 2 

^5 



(3.22) 

(3.23a) 
(3.23b) 



Note that there is one diffusive and one propagating 
sound mode in the incompressible limit when ^2 = 0. 
The attenuation (the imaginary part of q) of the diffu- 
sive modes is proportional to ^/uj and that of the sound 
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modes is proportional to to. This is the explanation of 
the large difference in attenuation of the two polariza- 
tions of transverse wave found by TW that makes 
NEs candidates for acoustic polarizers. We will return to 
this issue in Sec. IIV CI 

Equations (|3.21al) and H3.23a|) show that one misses the 
slow diffusive modes if one neglects the Frank energy. It 
is a legitimate question, however, under what conditions 
these modes can be observed experimentally. Manifestly, 
they should be observable directly in the symmetry di- 
rections in which the respective sound velocities vanish 
whereas they should not be seen in directions which differ 
significantly from these symmetry directions. Between 
the two extremes there will be a crossover from slow dif- 
fusive to propagating behavior at certain crossover an- 
gles. These angles can be estimated by comparing the 
magnitude of the terms in Eqs. H3.19() and (|3.21a(l or re- 
spectively (I3.23al) . For « we get for example from 
Eqs. H3.19|l and (|3.21a|l that the crossover is expected at 
an angle such that 




(3.24) 



The wave vectors in light scattering experiments on NE 
dynamics typically have a magnitude q ^ 10^ cm. The 
elastic moduli Ci and C4 should be comparable to the 
shear modulus of rubber, Ci ~ C4 ~ 10^ dynes/cm^. 
Hence, we estimate 

|sin6lo| ~ 10"*, (3.25) 

i.e., the range around the nematic direction in which the 
slow diffusive behavior is observable is extremely nar- 
row. Applying the same reasoning to the slow diffusive 
T-mode for w 7r/2 we obtain a further crossover angle 
6*^/2 with 

|cos0,/2| «10-^ (3.26) 

signaling the same extremely narrow angle range for slow 
diffusive behavior as above. 



IV. DYNAMICS WITH DISPLACEMENTS AND 
DIRECTOR 

Nematic elastomers, like their conventional nematic- 
liquid- crystal counterparts, are characterized by a Frank 
director n that responds dynamically to external forces. 
In this section, we will use the Poisson-bracket approach 
to derive phenomenological equations for the dynamics 
of both the director and displacements in nematic elas- 
tomers. Though our derivation is different from that 
of TW, our dynamical equations are in fact identical to 
theirs if they include contributions from the Frank free 
energy for the director. Our equations also reduce to 
the standard equations of nematohydrodynamics |33l l3^ 
when elasticity due to network crosslinking is turned off. 



It must be emphasized, however, that our equations pre- 
dict non-hydroydynamic modes characterized by a decay 
time r that does not approach infinity with vanishing 
wavenumber q. In any real system, there are many non- 
hydrodynamic modes with characteristic decay times Tq, . 
For a dynamical theory to provide a correct description of 
a system over a frequency range from zero to some max- 
imum frequency lom, it must include contributions from 
all modes, both hydrodynamic and non-hydrodynamic, 
with characteristic frequencies up to a few times ojm- 
An isotropic rubber is characterized by rather lo ng d ecay 
times Tpi of Rouse-like modes of chain segments 14^ . At 
frequencies uj such that wtr <C 1, its viscosities become 
frequency independent described by the hydrodynamical 
equations of an isotropic solid. However, when loth > 1, 
viscosities develop nontrivial frequency dependence, and 
hydrodynamics breaks down. In our theory for nematic 
elastomers, we assume that there is a single director re- 
laxation time T and that it is much larger than so that 
we do not need to worry about the frequency dependence 
of viscosities arising from Rouse modes. We will, how- 
ever, point out where these assumptions, also made by 
TW, can be modified. If r < r^, dynamics at frequencies 
LUTii > 1 will be dominated by Rouse modes, and it may 
be difficult to distinguish director relaxation modes from 
Rouse modes in experiments in which displacements are 
probed. The hydrodynamic description of the preceding 
section, however, remains valid when coTfj <^ 1. 



A. Equations of motion 

When we keep the director as a dynamical variable, 
our formulation is closely related to nematodynamics. 
The equation of motion for the director has a reactive 
coupling to the velocity v = u arising from the Pois- 
son bracket of ni{x) with gj{x')- This Poisson bracket 
is generally derived in Eulerian coordinates [33. in 
which there is a contribution v • Vn to the equation for 
dtii. This term can be combined with dtn to yield the 
Lagrangian time derivative n. The remaining part of the 
Poisson bracket is 

{n,(x),5,(x')} = -kjkdkSi^^^'). (4.1) 

The properties of the tensor Xijk are dictated by three 
constraints: First, the magnitude of the director has to 
be conserved, i.e., n • h = implying riiXijk = 0; sec- 
ond, the equations of motion must be invariant under 
n ^ — n implying Xijk must change sign with n; and 
third, under rigid uniform rotations, the director has to 
obey n = i(V x u) x n. The only tensors and vectors 
available to use for the construction of Xijk are Sij, rii, 
and the Levi-Civita anti-symmetric tensor Cijfc, and the 
only tensors that satisfy the first two conditions above 
are Sf^Uk and Sfj^nj where Sfj = Sij — ninj. Thus, the 
first two conditions imply that Xijk has two independent 
components, which can be expressed as parts symmetric 
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and antisymmetric under interchange of j and k: 



X 



X 



ijk 



Sik nk) 



A2 
2 



{Sf^n,~Slnk) ,(4.2) 



The third condition imphes that the coefficient A2 of the 
antisymmetric part must be equal to —1. There are no 
constraints on the value of A, which is equal to the ratio 
of two dissipative coefficients of the Leslie-Eriksen the- 
ory ■ When its absolute magnitude is positive, it de- 
termines the equilibrium tilt angle of the director under 
uniform shear [27l |. 

Using the Poisson brackets for the director and those 
for the momentum density and displacement discussed in 
the preceding section, we obtain the equations of motion 
for a nematic elastomer 



rii = Xijk OjUk - 1 7 — , 



5n, 



1 

P' 



■9i 



gi = Xkj'i. o-j- 1- Vijki OjdiUk . 



(4.3a) 
(4.3b) 
(4.3c) 



The first terms on the respective right hand sides of 
Eqs. H4.3a|l and (|4.3b|l as well as the first and second 
term on the right hand side of Eq. H4.3c|l stem from the 
Poisson brackets. H is the full elastic energy of NEs as 
stated in Eq. H3.1|l . The second term on the right hand 
side of Eq. H4.3a|l is a dissipative term that describes dif- 
fusive relaxation. There is no dissipative contribution to 
because NEs are tethered. 



Eq. Ij4.3b|l because NEs are tethered. Vijki is a viscosity 
tensor that has the same structure as rjijki- It has five 
independent components vi to 1^5. We use different sym- 
bols here for the viscosities than we used in Sec. IIIII so 
that we can cleanly keep track of differences between the 
two theories. 

We have assumed in Eqs. (|4.3(l that both F and i^ijki are 
local in time or, equivalently, that their temporal Fourier 
transforms are independent of frequency as they are at 
frequencies uj less than their respective inverse character- 
istic times TjT^ and . At larger frequencies, however, 
both F and Vijki do depend on frequency. In polymer gels, 
viscosities are proportional to ^/uJ at frequencies ^ 1 
because of the large number of closely spaced modes con- 
tributing at these frequencies. Though there are to our 
knowledge no microscopic calculations of F(ix'), there is 
no reason why many closely spaced modes should not lead 
to Rouse- like behavior at frequencies wTp ^ 1. There is 
also no reason why rr and tr should not be equal or 
nearly so. If the decay time t for the nonhydrodynamic 
director modes predicted by the the theory with the low 
frequency approximations to F and Vijki is greater than 
Tn and Tr, then this theory provides a correct description 
of the dynamics for frequencies < a; < min(r^^, TjT^) 
including uj . If on the other hand r < min(Tfl, rr), 

the full frequency dependence of F must be used for fre- 
quencies Lj > min(T^^, TjT^). In all cases, however, the 
hydrodynamic theory of the preceding section is regained 



in the limit uj ^ mm{T~^ , t]^^ , t^^) ■ To keep our dis- 
cussion simple, we will continue to assume that both F 
and Vijki are frequency independent. Our equations can, 
however, incorporate the frequency dependence of these 
quantities merely by replacing them by their frequency- 
dependent forms, F(a;) and 

At this point we would like to comment on the precise 
relation of Eqs. H4.3|l to the equations of nematodynam- 
ics. One retrieves the latter from the former by setting all 
the elastic constants in H, except for the Frank constants 
Ki, K2, and K3, equal to zero. Thus, we are guaranteed 
to obtain the well known modes of nematic liquid crystals 
in this limit. One can also make the observation that in 
this limit Eqs. 1)4. 3|) depend just on v = u rather than on 
V and u. If one seeks modes in terms of u rather than v 
one finds extra zero frequency modes that are spurious. 

Next we condense the equations of motion into a re- 
duced set of effective equations. We can simplify the 
subsequent analysis at the onset by eliminating either u 
or g with help of Eq. I)4.3b|l . To facilitate in contact 
with the work of TW we choose to keep u. For the same 
reason we opt to work with Q rather than n. This rep- 
resents no difficulty since Q and n are simply related via 
Eq. I|3.8|l . Collecting we obtain after some algebra the 
following effective equations of motion for Q and u: 



{[dt + rDi]6ab-rMab{'^)}Qb 



(4.4a) 



and 



pdtUa = --^^ dtd^lQa - XUaz\ + [C2 -|- V2dt\daU:, 



+ [C3 -I- Vzdt\daUii -I- 2[C4 + V^dt\dbUab 

. D2-DI 



C5 + Vbdt ~ ^ 



pdtUz 



X - 1 
2F 



dtdb\ 



dz 



Xuh 



dzQa , 

(4.4b) 



-I- [Ci -I- i^idt + C2 + v2dt]dzUz 
+ [C2 + i^2dt + C3 + vs.dt]dzUi, 



C5 + v^dt + ^ 



dbUbz 



D1+D2 



dbQb ■ 
(4.4c) 



To make contact to uniaxial solids we now take a brief 
detour and consider the simplified case A/a6(V) = 0. In 
this case Eq. (I4.4a|) is diagonal in frequency space and 
readily solved with the result 



VjJT2 



Di 1 — iuJTi 



(4.5) 



In writing Eq. H4.5|l we have used the relaxation times 
Ti = l/{TDi) andT2 = -A/(Fi?2). n « IQ-^s 
is essentially the relaxation time of the director. Pure 
hydrodynamic behavior is obtained when both luti and 
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UJT2 are much less than 1 as our calculations will ver- 
ify. In terms of the relaxation times, A is given by 
A ~ —T2D2/ (tiDi). Note that A is usually negative, 
but if it is positive there no inconsistencies arising from 
a negative T2. 

Inserting expression (|4.5(l into the time-wise Fourier 
transformed equations of motion for Ua and Uz , we obtain 



+ Ci{uj)dzUa 



(4.6a) 



-pu'^Uz = [Ci{uj) + C2{uj)]dzUzz + [C2(w) -f Ci{uj)]dzUi, 



+ Cfiu)dbU,z 



(4.6b) 



A2\ Dj (l-iwT2)2 



with Ci{lu) = Ci — iujiyi, 6*2(1^') = C2 — iujh'2 etc. and 
where 



2V ) 2Di l-iujTi 
- iojv^ + O {u^) , (4.7) 



2V 



1 - — 



(4.8) 



Note that Eqs. I|4.6|l are identical in form to Eq. (|3.15|l 
for uniaxial solids with K^^ — K§ — 0. Note also that we 
can identify 771 = vi and so on. Indicating the consistency 
of our two approaches, we can identify the renormalized 
viscosities rj^ and v^. As pointed our earlier, Rouse 
modes will be important for wr^, > 1 and wrr ^ 1 lead- 
ing to a frequency dependence of the viscosities. In this 
regime we have to let Vi —f Vi{u!) and T — > T{uj) . For 
ujTj^ 3> 1 and LOTr ^ 1, in particular, the viscosities are 
proportional to -^/tJ |44|. 

Before we move on we point out that Eq. (|4.4|) become 
identical to the equations of motion by TW if we take the 
incompressible limit and if we neglect the Frank energy, 
provided of course, that we take into account that some- 
what different conventions are used and provided that 
corresponding quantities are identified properly. A de- 
tailed comparison to TW at the level of final results for 
the modes will be given in Sec. IIV CI 



B. Mode structure 

1. The t-direction 

By applying Sba Qb (including the summation over a) to 
both sides of the Fourier transformed version of Eq. (|4.4a|) 
we obtain a diagonalized equation of motion for Qt that 
is readily solved with the result 



D2 1 - iu;T2 + q^Kt/D2 1 
Di l-iojTi + q^Kt/Di 2 



7: ^qzUt 



(4.9) 



where we have used the abbreviated notation 



q^Kt = K2ql + K^ql 



(4.10) 



Application of the same procedure to Eq. Ij4.4b|l yields 



A-hl 



-LO put 



■iLoqz 



+ 



2r 

D2-D1 . 



iQt 



A 



■ qzUt 



iqzQt 



C4{L^)q'iut 
D2 



qlut . 



(4.11) 



Ut is a purely transverse displacement that does not cou- 
ple to any displacement in the plane containing q. As a 
result, the elastic constants C2 and C3, which couple to 
u;, do not appear in this equation of motion. 

Next we insert Eq. (|4.9I) into (|4.11() which gives an ef- 
fective equation of motion for ut alone. This equation of 
motion has solutions with 7^ for frequencies satisfy- 
ing the secular equation 



uj'^p — Ci{ijj) q\ + 



A(A + 1) - Xt^^ ~ 2TC^{uj) 

4r 



Di [1 - iujTi + q^Kt/Di] 

A + 1 — r, ^ — Ato ^ o , T , -, 

1 ■^qlD2[l- iuJT2 + q^Kt/D2] = . 



4r 



(4.12) 



Evidently this secular equation is too complicated to find 
useful closed solutions for u. However, since we are in- 
terested in the long- wavelength behavior, we can pertur- 
batively determine solutions in the form of a power series 
in the wavevector. 

Let us first consider the simplified case of vanishing 
C's and D's. In this case we obtain as expected the well 
known slow and fast t-direction modes of nematic liquid 
crystals 



ujt,s = -I [K2 ql + i^ag'] \ r 



{l + >^?ql 
2{v^ql + 2viql)\ ' 



2vi q\ + 1^5 ql 
2p 



(4.13a) 
(4.13b) 



as well as the previously announced spurious zero- 
frequency mode. In writing Eqs. H4.13|l we have con- 
sidered as usual the hmit KpT/u <^ 1, where K stands 
symbolically for all the Frank constants and u stands am- 
biguously for all the viscosities. In comparing Eqs. (|4.13() 
and H4.19|l to the original results on nematic liquid crys- 
tals as given in Ref. jsj one should be aware of slight dif- 
ferences in the notations. If we wish to use the notation 
of Ref. im we have to replace in Eqs. (|4.13() F — > 7"^, 
vi 2vi, 1^4 V2, 1^5 2i^3, q± qi and q^. 

Now to the full secular equation H4.12|l . Solving 
this equation perturbatively leads also to three modes, 
namely one massive mode and two propagating modes. 
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The massive mode has the frequency 



-IT, 



2p 



(4.14) 



Note that Eq. (|4.14|) does not depend on C^. Hence, 
this mode is shared by soft and semi-soft NEs. As long 
as the sound velocities of the soft modes are finite, their 
frequencies are 



LOt,± = ± 



/ 2C4 



2p 



I - 



(4.15a) 



in full agreement with Eq. (|3.19a|l when Crf = and the 
identifications 1^4 = 774 and v^" = rj^ are made. For soft 
NEs the soft i- modes become diffusive ii q± = 0. The 
frequencies of these diffusive modes are easily identified 
with those given in Eqs. H3.20(l and H3.21|l . 

2. The T-direcUon 

Except for the i-direction an analysis of the modes is 
prohibitively complicated unless one resorts to the in- 
compressible limit. In this limit = 00 and ui — 0, 
and as a result, C2, C3, 1^2 and 1/3 do not appear in the 
equations of motion for ut- 

We need the equations of motion for the _L direction as 
intermediate results for studying the T-direction. Apply- 
ing qa (including the summation over a) to the Fourier 
transformed counterparts of Eqs. H4.4a|l and (|4.4b|l we 
obtain 

q^K^/D2 1 . 



Q± = - 



D2 1 



IUJT2 



Di 1 - iujTi + q^K^/Di 2 

D2 1 - iuJT2 - q^K^/D2 1 



-LU pU± 



Di 1 — iujTi 
A + 1 . 

2C4,{Lu)qlu 



q^K^/Di 2 
2 (9^^;^ 



iqzU± , (4.16a) 



qzU±) 



D2-D1 
2 

qziq^Uz 



iqzQj 



qzUA 



where we have used the shorthand 

Kiql+K^ql 



q^K^ 



(4.16b) 



(4.17) 



In addition to Eqs. H4.16|) we also need the Fourier trans- 
formed version of Eq. H4.4c() in the incompressible limit. 



'UJ^pUz 



A- 1 

2r 



■ iujqj_ 



2 (9^"^ 



- qzUi_) 



Ci (uj) qluz 
C5(c^) - 



D1+D2 



£1 

2 



2 

qj_{qj_u 



iq±QA 



+ qzui_). (4.18) 



Next we set up an equation of motion for ut that 
depends on u±, Uz and Qa_. Then we eliminate Q±^ 
with help of Eq. H4.16a|l . Finally, we exploit that = 
—qzUT/q and Uz = qiJu-Tlq in the incompressible limit. 
These steps provide us with an effective equation of mo- 
tion for Ut alone. In order to allow for solutions ut 7^ 
the frequencies have to satisfy a condition analogous to 
Eq. (|4.12|l . We opt not to write down this secular equa- 
tion here because it is rather lengthy and because it can 
be obtained in a straightforward manner from the ingre- 
dients given above. 

We proceed as above and first consider the simplified 
case of vanishing C's and Z?'s. As anticipated we obtain 
a spurious zero-frequency mode as well as the slow and 
fast T-direction modes of nematodynamics, 

UT.s = -i [Ki q\ + Kj,qf\ 

[q^-X{ql-ql)f 



X F 



2v^{ql~qlY + A{u^ + 2ui)qlqll' 

(4.19a) 

2{vi^2vi)qlql + V5 {ql-qlf 
^Tj = -i 77—2 . (4-19b) 

^ p q 

where KpT/v ^ 1 is implied. 

By perturbatively solving the full secular equation we 
extract the three T-direction modes for NEs. We find 
one massive mode 

WT,™ = -Tf' - iT [Ki ql + K3] ql 



2p 



(4.20) 



and two soft modes. For non- vanishing sound velocities 
the soft modes have frequencies 



/2(Ci+2C4)gi<z2 + cf (gi-g2 



2pq^ 



— I ■ 



. 2{ui + 2vi) qiqi + vi^ {qi - qj 



Apq^ 



(4.21) 



When = this result reduces, provided that the vis- 
cosities are properly identified, to Eq. (|3.19bp . In the 
case of ideal soft elasticity the sound velocities vanish if 
= or = 0. For g_L = the frequencies of the 
then diffusive T-modes are identical to those for the dif- 
fusive modes, see Eqs. (|3.20|) and (|3.21|) . For = we 
retrieve Eqs. and 



C. Comparison to TW 

Now we compare our findings to those by TW. First, 
we will demonstrate that our equations of motion H4.4(l 
are identical to the equations of motion by TW if we 
restrict ourselves to the incompressible limit and if we 
neglect the Frank energy. Second, we will compare our 
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TABLE L Correspondence between quantities used by TW 
and quantities used in our work. 



TW 


this paper 


TW 


this paper 


2Ci 


Ci 


2 A4 




2C2 


C2 


4 As 


i/5 + AV(2r) 


2C3 


C3 


71 


1/r 


2C4 


C4 


72 


-A/r 


4C5 


C4 


LU 




2 Ax 






Q2 


2A2 








2 A3 


«^3 







final results for the modes to those found by TW. Of 
course, we must take into account differences in conven- 
tions and notations in these comparisons. For guidance 
in identifying corresponding quantities, see Tabled 

In order to compare our equations of motion (|4.4(l to 
the equations of motion by TW we rewrite Eqs. (|4.4b() 
and H4.4c|l as 



-pdtu 



(4.22) 



where the aij are the components of the stress tensor 
g_. The specifics of the cry are easily gathered from 
Eqs. (j4.4b|) and (|4.4c|l . 



Cafc 



' zz 



2 [C4 + V^di\ Uab , 

[Ci + vidt] u^^ , 

2 ' 
A + 1 



C5 + i^sdt - - 
D2 + D1 



A(A + 1) 



2r 



dt 



r 

2 

A- 1 



2r 



r 



■dt 



(4.23a) 
(4.23b) 



(4.23c) 



(4.23d) 



where we have taken the incompressible limit. Equa- 
tions H4.23|) show clearly that g_ is not symmetric and 
that its antisymmetric part is 



(Taz - O-^a = i { [t^ ^ + dt] Qa - A[r2 ^ + dt] Uaz } 



(4.24) 



where we have used Eq. (|4.4a|) to obtain the last equality. 
Equation (|4.24ll reveals that the stress tensor try defined 
in Eqs. (|4.23|) is symmetric only when the Frank energy 
can be ignored. Of course, when the Frank energy is 
includedjh is al way s possible, following the procedures 
in Refs. |33| and [33|, to define a symmetric stress tensor 
afj yielding the same equations of motion as Uy . 

If we neglect the Frank energy, our stress tensor be- 
comes identical to the symmetric stress tensor used by 
TW provided that we take into account Tabled More- 
over, as easily can be checked, Eq. H4.4a|l becomes iden- 
tical to the balance of torques equation under these cir- 
cumstances. Therefore, our equations of motion (|4.4|l are 



equivalent to the equations of motion by TW provided 
that we take the incompressible limit and neglect the 
Frank energy. 

At this point it is interesting to compare the stability 
conditions which arc implicit in the dissipation function. 
Our equations imply in the incompressible limit that vi > 
0, 1^4 > and 1^5 > 0. While the first two conditions can 
be readily identified with the conditions ^1 > 0, A4 > 
in the work of TW, the situation is less obvious for the 
last condition. Note that we can re-express > in 
terms of the quantities used by TW as 



1^5 



7? 

271 - 



(4.25) 



which has an identical counterpart in the work of TW. 
Stated in terms of the relaxation times, this condition 
requires that 



(4.26) 



for ideally soft NEs. 

Next we come to the comparison of the results for the 
modes. Our first observation is here that the massive 
modes tot.m and LUx.m in unaccounted for by TW. Note 
that the soft modes in the t- and in the T-direction are 
referred to in TW as qSH waves and qSV waves, respec- 
tively. 

Let us first compare the findings regarding the t- 
direction. To foster this comparison we recast our re- 
sult (|4.15a|l as 



OJ 



1 



(4.27) 



To the order we are working, i.e. to order 0{q'^), the 
solutions to Eq. (|4.27l) and our ujt,m coincide. Note that 
Eq. H4.27|l is essentially identical to the t-direction secular 
equation for uniaxial solids. This can easily be checked 
by starting with Eq. (|4.t)a|) and by then switching to the 
t-direction. 

To make further contact to WT we eliminate the vis- 
cosities in Eq. H4.27|l in favor of the relaxation time 
TR « i^i/Ci « ••• « jy^/Cs « K + AV(2r)]/C5. The 
subscript R indicates that tr w 10^^ — lO^^s is of the 
order of the Rouse time of the polymers constituting the 
rubbery matrix. We obtain 



1 



2Ci ql + Ci[Lo) q^^ (1 - iuTR) , (4.28) 



where 



1 - lUTR 



2Dl (1 — ILUTl) (1 — ILOTr) 



(4.29) 



Taking into account Table we see that our C§{lo) is 
identical with the renormalized form of C5 found by TW 
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and that Eq. (|4.29() is in full agreement with the disper- 
sion relation of TW for the qSH waves. 

Now to the T-direction. To order 0{q^) our results for 
ut ± coincide with the solutions of 



L0^ = 



1 

2^ 



[C1H + 2C4H] 



9 9 



91 



(4.30) 



Note that Eq. H4.30|l is essentially identical to the 
T-direction secular equation for conventional uniaxial 
solids. TW considered the T-modes only for — 0. 
For Qz = Eq. 14.30|l reduces to 



UJ 



2p 



C^{cu)ql{l~tUJTR). 



(4.31) 



Using Table U find that Eq. (|OT|) is in full agreement 
with the result of TW for the qSV waves. 

Before summarizing our findings, we now briefly return 
to the property of NEs that makes them, as pointed out 
by TW, candidates for acoustic polarizers, viz. the large 
difference in attenuation between the t- the T-modes in 
the symmetry direction where qz — 0. For the <-direction 
it follows immediately from Eq. H4.28|l that the attenua- 
tion is proportional to co if Crf = and qz — 0. Equa- 
tion H4.31|l has two solutions for — 0, namely u; ~ Q 
and, with the proper identification — i]^ , the fast dif- 
fusive mode Eq. Ij3.23b|l . If the Frank energy is taken into 
account, the lo — mode becomes the slow diffusive mode 
of Eq. I|3.23a() . Due to their diffusiveness, the attenua- 
tion of these T-modes is proportional to ^/uj and hence 
much larger for low frequencies than the attenuation of 
the propagating t-mode. This large difference in attenu- 
ation can be used, in principle, to split the T-modes from 
the t-modes. 



fective shear modulus for shears in planes contain- 
ing the anisotropy axis respectively vanishes or is very 
small. In this paper, we have explored the dynamical 
consequences of this elasticity. We derived dynamical 
equations, involving only the displacement, valid in the 
hydrodynamic limit in which frequencies and wavenum- 
bers are respectively small compared to all characteris- 
tic microscopic inverse times and lengths in the system, 
and we determined that their mode structure is identi- 
cal to that of columnar liquid crystals in the soft limit 
when — 0. We also used the Poisson-bracket ap- 
proach to derive dynamical equations, which contain non- 
hydrodynamic modes, for both director and displacement 
and verified that they reduced to those derived by Ter- 
entjev, Warner, and coworkers [lll22j| when contributions 
from the Frank free energy can be ignored. We analyzed 
the mode structure of these equations assuming a sin- 
gle relaxation time for the director, which we took to 
be longer than any other characteristic decay time such 
as the Rouse time tj^. Our equations, however, permit 
the introduction of frequency-dependent dissipative co- 
efficients valid at frequencies higher than these inverse 
decay times. 

Rheological experiments at zero wavenumber have re- 
ported frequency-dependent storage and loss moduli that 
are in agreement with the predictions of the semi-soft 
theory that includes the director It would be in- 

teresting to map out the modes of nematic elastomers 
directly using light scattering. It may, however, be dif- 
ficult to access the true hydrodynamic limit because it 
applies in current NEs only for frequencies of order 100 
Hz or less and because inherent sample inhomogeneities 
may lead to extra scattering that could mask the signals 
of the characteristic modes of a homogeneous system. 



Acknowledgments 



V. SUMMARY 

Nematic elastomers exhibit the remarkable phe- 
nomenon of soft or semi-soft elasticity in which the ef- 



We gratefully acknowledge support by the Emmy 
Noether-Programm of the Deutsche Forschungsgemein- 
schaft (OS) and the National Science Foundation under 
grant DMR 00-95631 (TCL). 



[1] W. Warner and E.M. Terentjev, Liquid Crystal Elas- 
tomers (Clarendon Press, Oxford, 2003) 

[2] H. Finkelmann, H. J. Koch, and G. Rehage, Makromol. 
Chem. Rapid Commun. 2, 317 (1981). 

[3] P. G. de Gennes in Liquid Crystals of One and Two- 
Dimensional Order, p. 231, edited by W. Helfrich and G. 
Heppke (Springer, New York, 1980). 

[4] E. M. Terentjev, J. Phys. Cond. Mat. 11, R239 (1999). 

[5] P. G. de Gennes and J. Prost, The Physics of Liquid 
Crystals (Clarendon Press, Oxford, 1993). 

[6] S. Chandrasekhar, Liquid Crystals (Cambridge Univer- 
sity Press, Cambridge, 1992). 

[7] J. Kiipfer and H. Finkelmann, Macromol. Chem. Phys. 



195, 1353 (1994). 
[8] H. Finkelmann and H. Wermter, Polym. Mater. Sci. Eng. 

82, 319 (2000). 
[9] H. Finkelmann, E. Nishikawa, G. G. Pereira, and M. 

Warner, Phys. Rev. Lett. 87, 015501 (2001). 
[10] L. Golubovic and T. C. Lubensky, Phys. Rev. Lett. 63, 

1082 (1989). 

[11] H. Finkelmann, I. Kundler, E. M. Terentjev, and M. 

Warner, J. Phys. II (France) 7, 1059 (1997). 
[12] G. C. Verwey, M. Warner, and E. M. Terentjev, J. Phys. 

II (France) 6, 1273-1290 (1996). 
[13] M. Warner, J. Mech. Phys. Solids 47, 1355 (1999). 
[14] O. Stenull and T. C. Lubensky, Eurohys. Lett. 61, 776 



13 



(2003). See also, X. Xing and L. Radzihovsky, Eurohys. 

Lett. 61, 769 (2003). 
[15] O. StenuU and T. C. Lubensky, to appear in Phys. Rev. 

E (eprint: cond-mat/0307686 ) . 
[16] X. Xing and L. Radzihovsky, Phys. Rev. Lett. 90, 168301 

(2003). 

[17] J. L. GaUani, L. Hilliou, P. Martinoty, and P. Keller, 

Phys. Rev. Lett. 72, 2109 (1994); J. L. GaUani, L. Hilliou, 

P. Martinoty, F. Doublet, and M. Mauzac, J. Phys. II 

(France) 6, 443 (1996). 
[18] S. M. Clarke, A. R. Tajbakhsh, E. M. Terentjev and M. 

Warner, Phys. Rev. Lett. 86, 4044 (2001); S. M. Clarke, 

A. Hotta, A. R. Tajbakhsh, E. M. Terentjev Phys. Rev. 

E 65, 021804 (2002). 
[19] S. M. Clarke, A. R. Tajbakhsh, E. M. Terentjev, C. 

Remillat, G. R. Tomlinson and J. R. House, J. Appl. 

Phys. 89, 6538 (2001). 
[20] P. Stein, N. Assfalg, H. Finkelmann, and P. Martinoty, 

Eur. Phys. J. E 4, 255 (2001). 
[21] A. Hotta and E. M. Terentjev, Eur. Phys. J. E 10, 291 

(2003). 

[22] E. M. Terentjev, I. V. Kamotski, D. D. Zakharov, and L. 

J. Fradkin, Phys. Rev. E 66, 052701(R) (2002); see also 

L. J. Fradkin, I. V. Kamotski, E. M. Terentjev and D. D. 

Zakharov, Proc. R. Soc. Lond. A 459, 2627 (2003). 
[23] P. I. C. Teixeira and M. Warner, Phys. Rev. E 60, 603 

(1999). 

[24] E. M. Terentjev and M. Warner, Eur. Phys. J. E 4, 343 
(2001). 

[25] See for example, H. Goldstein, Classical Mechanics 
(Addison- Wesley, Reading, 1980). 

[26] J. L. Ericksen, Arch. Ration. Mech. Anal. 4, 231 (1960); 
Trans. Soc. Rheol. 5, 23 (1961). 

[27] F. M. Leslie, Q. J. Mech. Appl. Math. 19, 357 (1966); 
Arch. Ration. Mech. Anal. 28, 265 (1968). 

[28] Groupe d'Etudes des Cristaux Liquides (Orsay), J. 
Chem. Phys. 51 816 (1969). 

[29] P. C. Martin, O. Parodi, and P. S. Pershan, Phys. Rev. 
A 6, 2401 (1972). 

[30] K. Kawasaki, in Phase Transitions and Critical Phenom- 
ena, vol. 2, Eds. C. Domb amd M.S. Green (Academic, 
New York, 1976). 

[31] See, e. g., P. M. Chaikin and T. C. Lubensky, Princi- 
ples of Condensed Matter Physics, (Cambridge Univer- 
sity Press, Cambridge, 1995). 

[32] P.C. Hohenberg and B.I. Halperin, Rev. Mod. Phys. 49, 



435 (1977). 

[33] D. Forster, T. C. Lubensky, P. C. Martin, J. Swift, and 
P. S. Persham, Phys. Rev. Lett. 26, 1016 (1971). 

[34] Dieter Forster, Hydrodynamic Fluctuations, Broken sym- 
metry, and Correlation Functions (Addison Welsey, 
Reading, Mass., 1983); D. Forster, Phys. Rev. Lett. 32, 
1161 (1974). 

[35] See for example, J. Z. Xue, D. J. Pine, S. T. Milner, X. 
L. Wu, and P. M. Chaikin, Phys. Rev. A 46 6550 (1992). 

[36] See for example, H. Riskin, The Fokker-Planck Equation: 
Methods of Solutions and Applications (Springer Verlag, 
Berlin, 1984), and Ref. [s^. 

[37] Poisson brackets of coarse-grained variables require a 
proper interpretation. Perhaps most naturally, they can 
be regarded as coarse grained versions of the correspond- 
ing microscopic Poisson brackets. 

[38] H. K. Janssen, Z. Phys. B: Cond. Mat. 23, 377 (1976); 
R. Bausch, H. K. Janssen, and H. Wagner, ibtd. 24, 113 
(1976). 

[39] C. DeDominicis, J. Physique C 37, 247 (1976); C. DeDo- 
minicis and L. Peliti, Phys. Rev. B 18, 353 (1978). 

[40] H. K. Janssen, in: Dynamical Critical Phenomena and 
Related Topics, Lecture Notes in Physics, Vol. 104, ed. 
C. P. Enz (Springer, Heidelberg, 1979); H. K. Janssen, 
in: From Phase Transition to Chaos, ed. G. Gyorgyi, I. 
Kondor, L. Sasvari, T. Tel (World Scientific, Singapore, 
1992). 

[41] J. Schmidtke, W. Stille, and G. Strobl, Macromolecules 

33 2922 (2000). 
[42] M. Schonstein, W. Wille, and G. Strobl, Eur. Phys. J. E 

5, 511 (2001). 

[43] Beyond harmonic order, the left Cauchy strain ten- 
sor Vij = {diUj -\- djUi -\- dkUidkUj) /2 rather than the 
usual right Cauchy strain tensor Uij — {diUj + djUi — 
diUkdjUk) /2 couples to the director. See T. C. Luben- 
sky, R. Mukhopadhyay, L. Radzihovsky, and X. J. Xing, 
Phys. Rev. E 66, 011702 (2002). 

[44] See, e.g., M. Rubinstein and R. H. Colby, Polymer 
Physics (Oxford University Press, Oxford, 2003). 

[45] H. Stark and T. C. Lubensky, Phys. Rev. E 67, 061709 
(2003). 

[46] P. G. de Gennes, in Liquid Crystals of One- and Two- 
Dimensional Order, edited by W. Helfrich and G. Heppke 
(Springer, Berlin, 1980). 



